Feasibility of measuring the Shapiro time delay over 
meter-scale distances 



S Ballmer^ 2 3, s Marka^ and P Shawhan^ 

^Syracuse University, Syracuse, NY 13244, USA 

^TAMA - National Astronomical Observatory of Japan, Tokyo 181-8588, Japan 
^LjQO - California Institute of Technology, Pasadena, CA 91 125, USA 
CN "^Columbia University in the City of New York, New York, NY 10027, USA 

^University of Maryland, College Park, MD 20742, USA 

'""5 E-mail: sballmergligo . caltech . edu, smarka@phys.columbia.edu, 

pshawhan@umd.edu 

Abstract. The time delay of light as it passes by a massive object, first calculated by Shapiro 
Q-( in 1964, is a hallmark of the curvature of space-time. To date, all measurements of the Shapiro 

^ time delay have been made over solar-system distance scales. We show that the new generation 

bX) of kilometer-scale laser interferometers being constructed as gravitational wave detectors, in 

particular Advanced LIGO, will in principle be sensitive enough to measure variations in the 
Shapiro time delay produced by a suitably designed rotating object placed near the laser beam. 
We show that such an apparatus is feasible (though not easy) to construct, present an example 
design, and calculate the signal that would be detectable by Advanced LIGO. This offers the 
\^ first opportunity to measure space-time curvature effects on a laboratory distance scale. 
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1. Introduction 



The general theory of relativity asserts that the familiar force of gravity is fundamentally a 
manifestation of the geometry of space-time, which is given curvature by massive objects. 
Many measurements of observable quantities which are influenced by space-time curvature 
have confirmed the validity of the Einstein Equivalence Principle — the foundation on which 
all metric theories of gravity are based — and the apparent correctness of general relativity 
as the specific metric theory [IJ. Nevertheless, because of the importance of space-time and 
gravity to our understanding of the universe, opportunities for improved precision and new 
tests continue to be sought Q. 

The most stringent tests of general relativity have been obtained from measurements of 
the time delay of an electromagnetic wave as it passes by a massive object. Although this 
effect is a natural consequence of the equivalence principle and the curvature of space, it was 
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not called out by Einstein and was first described by Irwin 1. Shapiro in 1964 in the context 
of timing radar pulses reflected from Venus and Mars at superior conjunction [3]. The extra 
one-way time delay due to general relativity is 



where is the mass of the Sun, c is the speed of light, G is the gravitational constant, ri 
and r2 are the orbital radii of the Earth and the target planet, and d is the distance of closest 
approach of the radar beam to the center of the Sun (assumed here to be much smaller than ri 
and r2). For a radar reflection from Venus at the edge of the Sun, this amounts to a one-way 
time delay of 0.12 ms. 

The time delay may be different in alternative theories of gravity. In the Parametrized 
Post-Newtonian (PPN) formalism for arbitrary metric theories [|4l[T]|, the time delay is 



In the (1+7) factor, the 1 reflects the gravitational red- shift and is fixed for any metric theory 
of gravity, while the 7 measures the net curvature of space produced by the mass along the 
path of the light or radio waves [[U. For general relativity, 7 = 1, reproducing Shapiro's result. 
The possible deviation of 7 from unity has been limited by radar ranging to Venus, Mars, and 
spacecraft, as well as by timing the radio pulses from isolated and binary pulsars [El [6l |7l [H ; 
see [2| for a more recent summary of measurements. The most precise measurement was 
made using the Cassini spacecraft on its way to Saturn with a Doppler technique, yielding 
7— 1 = (2.1 ± 2.3) X 10~^ The same (1 + 7) factor governs the angular deflection of 
hght passing near the Sun, and very-long-baseline radio interferometry (VLB I) has been used 
to measure 7— 1 = (—1.7 ± 4.5) x 10""^ [fTOll . Lunar laser ranging also is sensitive to 7 in 
combination with other PPN parameters [fTTI . 

The Sun produces a measurable time delay at its edge despite the fact that most of its 
mass is located hundreds of thousands of kilometers away from the path of the radio beam. 
We may wonder whether smaller objects produce time delays for smaller separation distances 
according to the same model, i.e. the same effective value of 7. The PPN formalism does not 
allow for such a scale dependence, but there is no experimental evidence one way or the other 
for short distance scales. What time delay magnitude could we produce on a laboratory scale 
using a :iilO-ton mass, say? As a back-of-the-envelope calculation, utilizing (|2]) with M = 
10'^ kg, ri = r2 = 1 km, and d = 0.50 m, we estimate that 5t = 8.2 x 10^'^^ s. Remarkably, 
variations in that time delay due to changes in d can be detected using a suitable laser 
interferometer — for instance, one of the Advanced LIGO [12] gravitational- wave detectors — 
as we will discuss in section |2] below. 

If the distance from the mass to the beam changes from (ifar to d^car, the resulting change 
in time delay is 



Thus, we can spin a massive object near the interferometer beam to modulate the time delay 




(1) 




(2) 




(3) 



and then coherently integrate the signal over a long time period. Figure [T] is a sketch which 
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Figure 1. To produce a time- varying Shapiro time delay in one arm of a large interferometer, 
we envision placing several rotating masses next to the evacuated tube which encloses the 
laser beam. This figure shows the side-by-side beams of the two interferometers at the LIGO 
Hanford Observatory. The axis of rotation is parallel to the laser beams. Each section of the 
rotating mass system is supported by a bearing assembly in a separate vacuum enclosure (not 
shown); they rotate synchronously to produce a coherent effect on the laser beam. 



illustrates this basic concept. In section |3j we discuss the design considerations for the 
rotating mass and present a specific example design which is sufficient to make an interesting 
measurement of the time delay, as described in section]?} We briefly examine several possible 
systematic effects and implementation issues in sections |5] and |6] and then conclude with a 
discussion of how future interferometers could make better measurements possible. 

2. Laser Interferometers 

Direct detection of gravitational waves [13J is the primary goal of modern interferometric 
gravitational wave detectors such as LIGO (Ml, VIRGO IBl [H, GEO600 [fTTll, 
TAMA300 [fTSll and CLIO [19| which have successfully collected a few years of high- 
sensitivity data so far. Second-generation gravitational wave detectors now being designed 
and constructed — namely Advanced LIGO lfT2ll . Advanced Virgo [|20l and the Large-scale 
Cryogenic Gravitational-wave Telescope (LCGT) [1211 — will greatly extend the reach for 
gravitational wave signals. Each of these detectors is capable of measuring extremely small 
changes in the effective length difference of its two perpendicular "arms", as would be induced 
by a passing gravitational wave. However, they can also be viewed as unique "laboratory" 
instruments capable of measuring induced displacements at the ~ 10~^^ m level, equivalent 
to time delays of ~ 3 x 10~^^ s, for integration times of order one year. 

This feature opens up new possibilities for fundamental science. While gravitational 
wave searches target signals with various durations, and either accurately-known or unknown 
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Figure 2. Nominal expected noise levels of second-generation gravitational wave detectors as 
a function of frequency, adapted from numerical models provided courtesy of the LIGO ll22ll 
and Virgo EqIIII projects, and published by LCGT ED. The Advanced LIGO curve is for 
the "zero-detuning, high power" configuration; see ll22ll for a brief summary of alternative 
configurations and caveats. The LCGT curve is for the "broadband RSE" configuration. The 
lowest strain noise is obtained for frequencies between about 70 and 300 Hz. Also shown is 
a representative noise curve for a possible third generation detector design ll24l . expected to 
have an order of magnitude lower noise and to have access to lower frequencies. 

waveforms, precision experiments in which the light travel time in one of the interferometer 
arms is modulated at a given frequency in a well-controlled manner can enjoy all the benefit 
of sensitive lock- in experimental techniques and extended integration times. 

Advanced LIGO will consist of three 4 km interferometric detectors at the two existing 
LIGO sites: the Hanford Observatory in Washington State (home to two detectors) and the 
Livingston Observatory in Louisiana. (A proposed alternative plan, to install one of the three 
Advanced LIGO detectors at the Gingin site in Australia, leaving just one at Hanford, is under 
consideration at the time of this writing.) The detectors will be installed beginning in 201 1 and 
are expected to be fully operational around 2015, ultimately reaching an order of magnitude 
lower strain noise than the initial LIGO detectors. 

The Advanced LIGO detector design has a beam splitter and suspended mirrors (test 
masses) at the ends of its orthogonal arms, as well as additional partially-transmitting mirrors, 
to form a power- and signal-recycled Michelson interferometer with Fabry-Perot arm cavities. 
A gravitational wave induces a time-dependent metric strain on the detector which changes the 
relative lengths of the arms. While acquiring scientific data, feedback to the mirror positions 
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and to the laser frequency keeps the optical cavities near resonance, so that interference in the 
hght from the two 4 km long arms recombining at the beam splitter depends on the difference 
between the lengths of the two arms modulated by the signal to be measured. A photodiode 
senses the light, and a digitized signal is recorded. Then the data are calibrated and converted 
into a strain time series. The detectors have a sensitive frequency band extending from a few 
tens of Hz to a few kHz, as shown in figure |2} limited at low frequencies by seismic noise 
and radiation pressure noise and at high frequencies by laser shot noise. Third-generation 
detectors are envisioned to push the limits of ground-based interferometry and will likely be 
limited at low frequencies by Newtonian gravity-gradient noise ll25l . 

3. Rotating Mass Design 

Since laser interferometers like LIGO have essentially no sensitivity at DC, it is necessary 
to modulate the Shapiro delay by changing the effective distance between the beam and the 
relevant mass. The easiest way to achieve this is a rotating mass system. In this section 
we describe one possible design for a device with two "arms" extending from the rotation 
axis that would provide a Signal-to-Noise Ratio (SNR) of ~8 for one year of integration with 
Advanced LIGO. (There is no plan at this time to actually construct and operate such a device, 
but it could be considered as an addition to Advanced LIGO or another interferometer in the 
future.) 

As indicated by ([3]), the geometry of the rotating mass setup enters only logarithmically 
into the Shapiro delay, while the mass enters linearly. The closest effective distance of the 
mass to beam, dncarj is constrained by their finite cross-sectional sizes. As a consequence it is 
hard to reduce the required mass by optimizing the rotating mass geometry. For the example 
design we assume a mass of L5 x 10^ kg for each of the two arms. 

One of the limiting factors is the tensile strength of the material holding the rotating 
mass together. This of course favors low spin frequencies, and therefore requires a laser 
interferometer with good low-frequency sensitivity. In the case of Advanced LIGO we would 
like to have the primary harmonic at about 50 Hz, well separated from known instrumental 
lines in the detector noise spectrum. Thus, a baseline design with two arms must rotate at 
25 Hz. Note that this is still much slower than the inverse light travel time for all relevant 
distances; thus a static calculation of the Shapiro delay is still appropriate. A rotating mass 
system design with more arms would have a lower spin frequency, but in order to preserve the 
size of the variation of the Shapiro delay, one would have to increase the radius by the same 
factor. Since the whole rotating mass system needs to spin in vacuum and it must fit next to 
the beam tube we choose 1.5 m as the maximal practical rotating mass radius and stay with 
the conceptually simpler two-arm design and a 25 Hz spin frequency, even though one might 
win a factor of two in required tensile strength by going to 4 arms. 

With these parameters the rotating mass structure needs to support a centrifugal force 
of about 6 X 10^ N. A simple rectangular slab of steel, tungsten or titanium is not strong 
enough, and only a fraction of the mass would pass close to the laser beam. A better way 
to build a rotating mass system with the necessary strength is to concentrate the mass near 
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the ends of the arms, with a light-weight central support structure for suspension, and an 
outer shell of carbon fiber composite. Figure [3] shows a sketch of this rotating mass design, 
including the key dimensions. With a carbon fiber composite cross-sectional area of 1 
(5 cm thickness on both sides and 10 m cumulative rotating mass length along the axis of 
rotation) this requires a tensile strength of 600 MPa. A composite such as HexTow IM9 from 
HEXCEL corporation |j26l, for example, provides a factor of 5 safety margin on this. 

Another constraint on the geometry of the setup, specifically the length of the 
interferometer arms, comes from the direct gravitational coupling from the rotating mass to 
the test masses. A good place to put the rotating mass is near the middle of one interferometer 
arm. The test mass acceleration at the fundamental spin frequency /spin cancels due to the 
symmetry of the rotating mass, while the acceleration variation at 2/spin has amplitude 

15 GMr 2 (r + d) 2 
" = y (L/2)a 

where M is the mass of one arm, r is the rotating mass radius, d its minimal separation 
from the laser beam, and L is the interferometer arm length. For a 25 Hz spin frequency and 
including the acceleration of both test masses of the interferometer arm the length change at 
/ = 50 Hz has amplitude 

240 GMr 2 (r + c/) 2 

d^=^ PL^ 

'384 meter\6 2GM 

(6) 



^384 meter y 



c2 



Since the size of the Shapiro delay is on the order of the Schwarzschild radius 2GMj(? of 
one of the arm masses, divided by c, any laboratory- scale interferometer (L ~ 40 m or less) 
would require an extremely precise cancellation mechanism. For the L = 4 km arm length of 
LIGO this coupling is down by more than 6 orders of magnitude. 



4. Signal Analysis 

For a given rotating mass, such as the one described in the previous section, it is 
straightforward to numerically calculate the gravitational potential V at every point in its 
vicinity. A light ray which passes near the rotating mass experiences a total time delay which 
is simply given by 

5t = J{l + -f)Uds (7) 

where the integral is along the path of the light ray [1]. The light also experiences angular 
deflection during its passage, but for this system the deflection is only of order ~ 10^22 j-^jj 
The slightly increased length of the path due to this deflection is quadratic in the angle, and 
thus has a negligible effect on the time delay calculation. Therefore we may take the light path 
to be straight in the z direction (along the laser beam) and calculate the integrated time delay 
as a function of its transverse position when it passes by the rotating mass. Also, because the 
potential is additive, the spatial extent of the rotating mass in the z direction does not matter; 
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Figure 3. Cross-section of the example rotating mass system. All units are in cm. The support 
structure centered on the rotation axis is adequate to support the static weight of the steel 
weights, which represent most of the mass of the system. The carbon fiber epoxy composite 
outer envelope provides tensile strength when the system is spinning. The center of the laser 
beam is located 30 cm from the outer edge of the masses. 



we can project all of its mass density to the x-y plane. Figure |4] shows contours of constant 
time delay (isochrones) around our example rotating mass design. 

The time delay experienced by the light in the interferometer beam at any instant depends 
on the orientation of the rotating mass. As seen from the point of view of the rotating mass, 
the interferometer beam travels in a circle around the axis, as illustrated by the dashed lines 
in figure |4] for beams with closest-approach distances of 30 cm and 80 cm, representing the 
two interferometer beams at LIGO Hanford. The time delay at the center of the beam as a 
function of time is shown in figure |5] for one full rotation of the rotating mass. This is the 
signal which would appear in the interferometric data stream. Because the form of the signal 
is known accurately, it may be extracted from the data using lock-in techniques. 

The signal is not purely sinusoidal because the distances from the different parts of 
the rotating mass are comparable to the size of the rotating mass. Figure [6] shows the 
Fourier components of the signal (for the nominal spin frequency of 25 Hz) superimposed 
on the expected sensitivity curve for Advanced LIGO assuming an integration time of one 
year. Based on this comparison, the fundamental frequency component at 50 Hz should be 
detectable with amplitude SNRs of 7.2 and 4.4, respectively, in the closer and farther beams. 
The first harmonic (at 100 Hz) should also be detectable in the closer beam with SNR of 2.3. 
Comparing these different measurements would test the consistency of the results with the 
expected form of the signal, while combining them in quadrature yields a total SNR of 8.7. 



Feasibility of measuring the Shapiro time delay over meter-scale distances 8 




Horizontal distance [cm] 

Figure 4. Shapiro time delay of the laser beam as a function of the position around the rotating 
mass. Note that the contour pattern rotates with the mass. The contour labels are in units of 
100 yoctometers (10~^^ m). The small black circles illustrate the relative size and position 
of the laser beams at the time of the closest approach. The large dotted circles trace out the 
positions of the centers of the laser beams relative to the rotating mass for a full revolution for 
the specific configuration discussed. The paths of the dashed circles through the isochrones 
map out the curves shown in figure [5] 

5. Possible systematic uncertainties 

To interpret the measured signal as being due to the Shapiro time delay, we must have 
sufficiently accurate knowledge of the geometry of the apparatus and must be able to rule 
out other possible systematic effects which could affect and/or fake the Shapiro delay signal. 
However, given that we only expect an SNR of ~8 above the instrument noise, we only need 
to limit these systematic effects to a few percent of the expected signal amplitude. 

First, we must know the separation between the rotating masses and the interferometer 
beam to within a few percent, i.e. perhaps 1 cm. The rotating mass system parts can be 
machined and assembled with much tighter tolerances, and also measured precisely after 
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Figure 5. Shapiro time delay of tiie laser beam as a function of tiie time for one full rotation of 
the rotating mass system. The two ttaces correspond to laser beam positions 30 cm and 80 cm 
away from the outer edge of the rotating mass. 



assembly. When spun up, the rotating mass system would stretch by about 5 mm due to the 
extreme tension on the carbon fiber composite outer layer; this should be taken into account in 
the data analysis and may not be known precisely (due to uncertainty in the bulk modulus of 
the carbon fiber composite) but would be a small correction. Knowledge of the interferometer 
beam position within its vacuum tube would be a greater challenge, but can probably be 
determined to 1 cm or better. The Advanced LIGO alignment system should keep the position 
of the beam stable to much better than 1 mm over the course of a year-long data run. There 
would also be a horizontal gradient in the time delay over the area of the interferometer 
beam, but the intensity-weighted mean delay (over the Gaussian beam profile) would be well- 
determined and the gradient is not expected to cause any other problems. Additionally, the 
diameter of the laser beam's vacuum tube may have to be reduced locally to allow the required 
close approach of the rotating mass to the laser beam. 

It may be advantageous to minimize the length of beam tube devoted to the rotating 
mass system by placing synchronized rotating units on both sides of the beam tube. This also 
would reduce the time delay gradient over the area occupied by the laser beam. In the case 
of LIGO Hanford, rotating mass systems on both sides would produce equal time delays in 
both interferometers (but would only partly reduce the time delay gradient across each laser 
beam). 

We must also consider parasitic coupling of the rotating mass system's motion to the 
interferometric signal through mechanisms other than the time delay. As mentioned in 
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Figure 6. Fourier components of the Shapiro time delay signal. The diamonds and squares 
correspond to laser beam positions of 30 cm and 80 cm away from the outer edge of the 
rotating mass. Also shown is an example Advanced LIGO sensitivity curve, scaled for 1 year 
of integration time. 



section [3| direct gravitational coupling to the test masses would be negligible due to the 
length of the interferometer arms. Electromagnetic interference from the motor driving the 
rotating masses should be negligible due to being so far away from the sensing electronics. 
Additionally one can use a gear so that the motor runs at a different frequency. 

Another concern is local vibration of the interferometer beam's vacuum tube, directly 
excited by the time-varying gravitational field of the rotating mass system. This vibration 
would have a significant component at twice the spin frequency and could couple into the 
interferometric signal by modulating light scattered off the walls of the beam tube. 

Baffles placed at an appropriate distance can sufficiently reduce the effects of light 
scattering. These baffles must be mechanically isolated from vibrations coming from either 
the beam tube or the ground. 

The residual coupling from such baffles can also be estimated using (|6]), with an 
additional scattering factor. 

, / Eseat \ /384 meter \ 6 2GM 

'""^ " \-^) ^ 

Thus placing the baffles a few hundred meters away from the rotating mass should provide 
sufficient isolation. Secondary scattering should also be analyzed, but that goes beyond the 
scope of this paper. 

Such a massive rotating mass system, if even slightly unbalanced, would vibrate the 
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ground as it rotates. We should be able to tune the balance below the sensitivity level of 
seismometers at initial assembly. However, residual vibration might evolve due to aging or 
slow differential stretching of the rotating mass's carbon fiber composite layer when it is spun 
up to very high speeds. If the vibration is not sufficiently attenuated as it travels through 2 km 
of ground, it could shake the sensing optics and introduce a spurious signal. However, the 
vibrational signal would be generated primarily at the spin frequency and should be detectable 
with seismometers near the rotating mass systems and near the interferometer test masses. 
Thus a vibrational signal carried through the ground would largely be distinguishable from 
the time delay signal, which appears at twice the spin frequency and harmonics thereof. 

We also note that any residual false signal due to vibrational coupling is likely to appear 
in the data with a phase different from what is expected from the relativistic time delay signal; 
this provides a useful consistency test. 

6. Engineering considerations 

Actual construction and operation of a rotating mass system like the one described in this 
paper would present several engineering challenges, but we beheve that they could all be met 
with existing technologies. While detailed solutions are beyond the scope of this exploratory 
study, we can at least comment on the major issues. 

In the example rotating mass system design, a total mass of ~30 tonnes with a total 
length of 10 m must be supported reliably and spun with little friction. Rather than attempting 
to do that with a single assembly 10 m long, we imagine that it would be divided into several 
independent sections, each supported by multiple bearings to be able to tolerate the failure of 
a bearing. The sections would be rotated by independent motors. The rotation of each section 
would be monitored with a rotary encoder and fed back to the drive motor to maintain the 
desired spin frequency and phase. The rotation and drive force would also be monitored for 
signs of possible drive or bearing failure. 

The kinetic energy of each rotating section would be of order 100-200 MJ, comparable to 
the energy in commercially available flywheel units used for energy storage ['27]. Obviously, 
a structural failure must be avoided at all costs; the safety of the facility and personnel would 
be crucial considerations if this experiment is to be carried out. 

7. Discussion 

In this paper, we have explored the possibility of measuring the effects of space-time curvature 
over laboratory distance scales for the first time. We have found that it is technically feasible 
to make this measurement using the Advanced LIGO interferometers. However, it is far from 
easy. We had to resort to the use of high-tensile- strength composite materials to obtain a 
rotating mass system with sufficient mass and rotation speed to produce a measurable effect 
with an SNR of ~8 integrated over a year of data. This is sufficient to measure 7 with a 
precision of ~25%, which could distinguish between the cases 7 = 1 and 7 = 0, for instance. 
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Future interferometers could improve this type of measurement, especially if their low- 
frequency sensitivity is better [12811 so that the mass can be rotated at a lower frequency. 
In particular, a design study is currently underway for a future underground "Einstein 
Telescope" L29il with projected sensitivity shown by the lowest curve in figure |2j That 
detector could support a Shapiro time delay measurement at a spin frequency of ~15 Hz 
that would be an order of magnitude better than the Advanced LIGO example presented in 
this paper. Besides the planned gravitational-wave interferometers, a dedicated interferometer 
could perhaps be used to make this time delay measurement. Although we cannot hope to 
approach the precision of the radar ranging measurements with this technique, it is interesting 
to be able to test gravity in an entirely different distance regime. 
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